Phase diagram evolution at finite coupling in strong coupling lattice QCD 
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We investigate the chiral phase transition in the strong coupling lattice QCD at finite temperature and density 
with finite coupling effects. We adopt one species of staggered fermion, and develop an analytic formulation 
based on strong coupling and cluster expansions. We derive the effective potential as a function of two order 
parameters, the chiral condensate a and the quark number density p q , in a self-consistent treatment of the next- 
to-leading order (NLO) effective action terms. NLO contributions lead to modifications of quark mass, chemical 
potential and the quark wave function renormalization factor. While the ratio p c {T — 0)/T c (/x = 0) is too 
small in the strong coupling limit, it is found to increase as = 2N c /g 2 increases. The critical point is found 
to move in the lower T direction as f3 increases. Since the vector interaction induced by p q is shown to grow 
as f3, the present trend is consistent with the results in Nambu-Jona-Lasinio models. The interplay between two 
order parameters leads to the existence of partially chiral restored matter, where effective chemical potential is 
automatically adjusted to the quark excitation energy. 

PACS numbers: 11.15.Me, 12.38.Gc, ll.10.Wx, 25.75.Nq 



I. INTRODUCTION 

Exploring the chiral phase transition and its phase dia- 
gram in Quantum Chromodynamics (QCD) is one of the 
most challenging problems in quark hadron physics. The 
chiral phase transition may really happen in compact astro- 
physical phenomena such as the early universe and compact 
stars, and can be investigated in heavy-ion collision experi- 
ments. The large magnitude of the elliptic flow parameter 
observed in the relativistic heavy-ion collider (RHIC) exper- 
iments indicates the formation of strongly interacting quark- 
gluon plasma (H 0] at high temperature. The future ex- 
periments at FAIR and in low energy programs at RHIC are 
expected to provide new discoveries in the phase diagram in- 
vestigations. 

The most rigorous and reliable framework to investigate 
the QCD phase transition would be the lattice QCD Monte- 
Carlo (MC) simulations. In the high temperature and low 
density region, the lattice MC can provide the quantitative 
predictions, and the critical temperature is estimated to be 
T c ~ (160 - 190) MeV fill]. In comparison, the lattice 
MC simulations do not work well in the high baryon den- 
sity region because of the notorious sign problem of the Dirac 
determinant. Many ideas have been proposed to overcome 
this problem |0], for example, the Taylor expansion around 
M = Hi, analytic continuation 12, 0, canonical ensem- 
ble method [11], improved reweighting method lfl2ll . and the 
density of states method loll . It has become possible to ac- 
cess the relatively small density region V/T < 1.0 1101. One 
of the interesting objects is the critical end point (CEP) lfl4ll . 
Recent works indicate that CEP may locate in the region 
fi/T > l.Q 03 EE E3], while de Forcrand and Philipsen 
suggest that CEP might not exist lfl8ll . The larger chemical 



* miura@yukawa.kyoto-u.ac.jp 



potential region is still under debate, and alternative methods 
are necessary to reveal the whole structure of the phase dia- 
gram. 

In the phase diagram investigation, the strong coupling lat- 
tice QCD (SC-LQCD), the lattice QCD formulation based 
on the expansion of the action in the power series of the in- 
verse bare coupling squared (1/ g 2 ), is an interesting approach, 
since the sign problem can be weakened or avoided. The 
SC-LQCD was applied first to the confinement study in pure 
Yang-Mills theories I2TI1 . Wilson suggested the con- 

finement mechanism in an analytic study of the strong cou- 
pling limit (SCL) of lattice QCD j2ll|. Creutz showed that 
the (3 = 2N c /g 2 dependence of the lattice spacing a in the 
MC simulation smoothly connects the strong coupling behav- 
ior and the continuum spacetime scaling behavior [ 22] . By us- 
ing the character expansion technique, Miinster demonstrated 
that the pure Yang-Mills SC-LQCD with high order correc- 
tions could explain the above MC results II23I1 . The scaling 
behavior in MC simulations indicates that the confinement is 
actually realized in the continuum spacetime, and the success 
of SC-LQCD suggests that the scaling region would be ac- 
cessible in SC-LQCD within the conversion radius, which is 
shown to be finite in pure Yang-Mills theories M24I1 . We may 
expect that the scaling and convergent properties are also kept 
with fermions. Then it would be possible that the SC-LQCD 
could provide useful results on the phase diagram in the whole 
region of the T — fj, plane. 



The SC-LQCD with fermions has a long history of chiral 
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and many theoretical tools have been developed; the large di- 
mensional or 1/d (d = spatial or spacetime dimension) expan- 
sion Q30Q, the finite temperature treatments in the Polyakov 
gauge [34] and in the tem por al gauge 1 36], the finite quark 
chemical potential effect 13511 with the help of the lattice 
chemical potential ll62"tl . The analytic expression of the SCL 
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effective potential has been derived at finite T l34l l36ll or 
at finite fi 13511 . In 1990's, phase diagram studies met some 
successes ll37l l38l l39ll based on the SC-LQCD effective ac- 
tion II 3 61] . We also find several works on the Polyakov 
loop 16311 based on SC-LQCD, and the functional form of 
the effective potential in SC-LQCD has provided basic ingre- 
dients in the SC-LQCD related models |64], [H |6(J and the 
Nambu-Jona-Lasinio model with Polyakov loop (PNJL) j67ll . 

Based on these successes, the SC-LQCD is recently revis- 
ited and expected to provide an instructive guide to QCD un- 
der extreme conditions. The pure Yang-Mills SC-LQCD ll23tl 
is extended to finite T B68I1 . and the shear viscosity below 
the deconfinement transition temperature is also studied ll69tl . 
In color SU(2) QCD, the interplay between the diquark con- 
densate A and the chiral condensate a is investigated at fi- 
nite T and \i B4011 . The diquark effect is also investigated for 
color SU(3) at zero temperature lT4lTl . We find remarkable 
develojjments in the SU(-/V C = 3) phase diagram investiga- 
tions ll42ll43ll44ll5lll . where a "naive" structure of the phase 
diagram with first and second transition lines separated by a 
tri-critical point (TCP) is obtained in the strong coupling limit 
with zero quark masses. With non-zero quark masses, TCP 
becomes a critical end point (CEP), whose discovery is one of 
the physics goals in low-energy programs at RHIC. 

In order to discuss the chiral symmetry on the lattice, 
the SC-LQCD has been developed in several fermion for- 
malisms. We find some pioneering works based on the 
staggered JH HI, the Wilson JH Hi] and the naive JUt] 
fermions. The domain-wall ll70ll and the overlap ifTTll fermion 
provide modern formulation of the lattice chiral symmetry, 
and some SC-LQCD based investigations are found in [58, 59] 
(domain-wall) and ll60l l6lll (overlap). In the present work, 
we adopt one species (n/ = 1) of (unrooted) staggered 
fermion. Its simple realization of the chiral symmetry on the 
lattice i25l l72i l73ll is useful to develop analytic formulations. 
It has been theoretically suggested 11321 17411 and numerically 
established [75] that the unrooted staggered QCD is equiva- 
lent to the four flavor (Nf = 4) QCD with degenerate masses 
in the continuum limit. 

In this paper, we investigate the phase diagram evolution 
with finite coupling effects. We employ one species of (un- 
rooted) staggered fermion, and take account of the next-to- 
leading order (NLO, 0(1/ g 2 )) terms in the strong coupling 
expansion. We concentrate on the leading order of the large 
dimensional (1/d) expansion [30] for simplicity. The gluon 
field is evaluated in the Polyakov gauge 113411 with respect 
for the finite temperature T effects, and the finite density ef- 
fects are introduced via the quark chemical potential [i on the 
lattice [62]. In these setups, we derive an analytic expres- 
sion of the effective potential in the mean field approxima- 
tion. In particular, the following points are newly developed. 
Firstly we introduce the NLO effective action terms through 
the systematic cluster expansion. Secondly we evaluate the 
NLO effective action by using a recently developed extended 
Hubbard-Stratonovich (EHS) transformation [45, 46]. As a 
result, several auxiliary fields including the chiral condensate 
a are introduced on the same footing, and the NLO effects 
are self-consistently evaluated. In particular, we find that the 



quark number density naturally appears as an order parameter, 
whose self-consistent solution in equilibrium plays essential 
roles in the large /i region. This point would be an advan- 
tage to the previous works with NLO effects US S, [H . 
Thirdly, we discuss the evolution of the first and second or- 
der transition lines and the critical point with (3 — 2N c /g 2 . 
The finite coupling effects on the critical point have not been 
investigated before. Fourthly, the NLO contribution is ex- 
pressed as modifications of the constituent quark mass, chemi- 
cal potential and the quark wave function renormalization fac- 
tor. Hence the mechanism of the phase diagram evolution be- 
comes clear. 

While we are working based on the strong coupling expan- 
sion with NLO effects, we expect the present work would give 
a valuable picture in understanding the QCD phase diagram 
in the real world through the relation with the MC simula- 
tions. Since the SC-LQCD is based on the same formula- 
tion as lattice MC simulations, its results should be consis- 
tent with MC results as long as the applied approximations 
are valid. This speculation is supported by previous works 
on the hadron mass spectrum J27l[30l l3lll . Very recently, the 
structure of the phase diagram suggested in SC-LQCD is qual- 
itatively confirmed by a lattice MC simulation in the strong 
coupling limit 115 ill based on a Monomer-Dimer-Polymer for- 
mulation ll53ll . Thus the phase diagram in the strong coupling 
limit is established from both side of analytic and numeri- 
cal studies, and provide a good starting point to explore the 
true phase diagram by evaluating finite coupling effects in the 
strong coupling expansions. We find MC studies using one 
species of unrooted staggered quarks, and the results around 
(3 ~ 5 have been extensively discussed ifiol, l76tl . In order to 
compare the SC-LQCD results with those in MC simulations, 
we discuss the results in the region (3 < 6 expecting that these 
(3 values are within the conversion radius. 

Although the number of flavors (Nf — 4n/ = 4) used in 
the present work is different from the real world (Nf = 2 + 1), 
we could provide valuable results for the phase diagram inves- 
tigations. Flavor dependence of the phase diagrams at strong 
coupling has been studied by using several species (nf = 2,3) 
of staggered fermions lf38l l43ll . and we find that the phase 
diagrams with n/ = 2 and 3 are qualitatively the same as 
that with rif = 1. The critical chemical potentials at T = 
(Mc,t=o) are found to be almost the same. The critical temper- 
atures at fx — (T CjM= o) are found to be around 1.2 ll38il43ll 
and 1.06 ll38tl for nf = 2 and 3, respectively. These values 
differ from the result of n } = 1 (T C:AI=0 = 5/3) by 30-40 %, 
but the obtained phase diagram structure is very similar. The 
flavor dependence of the phase diagram is found to be mod- 
erate also in the continuum region. In Ref. [10], the small /1 
region of phase diagram is investigated by using MC simu- 
lations with four flavor staggered quarks, and the results are 
compared with two flavor results 19J. The difference between 
the phase boundaries in two and four flavor cases is at most 7 
% within a region N c fi < 500 MeV [10]. Thus, a "shape" of 
the phase boundary may not be crucially affected by the flavor 
effects. It should be noted that the number of flavors is impor- 
tant to some of the key features of the phase diagram, such as 
the order of the phase transition and the position of the critical 
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point. 

The organization of this paper is as follows. In Sec. [II] we 
provide a brief review on the strong coupling (l/g 2 ), the large 
dimensional (1/d) and cluster expansions, and derive the ef- 
fective action including the 0(l/g 2 , l/d°) effects. In SecHlH 
we derive an analytic expression of the effective potential. 
In Sec. [IV] we investigate the phase diagram evolution with 
(3, and focus on the mechanisms of the critical temperature 
and chemical potential modifications. The "partially chiral re- 
stored (PCR) matter" is found to appear in the high density 
region, and we also discuss its origin. Finally we summarize 
our work in Sec.[V] All through this paper, we use the lattice 
unit a = 1, and physical values are normalized by a. 



discrete values, we consider T as a continuous valued tem- 
perature. We take account of finite T effects by imposing 
periodic and anti-periodic boundary conditions on link vari- 
ables and quarks, respectively. We take the static and diago- 
nalized gauge (called Polyakov gauge) for temporal link vari- 
ables with respect for the periodicity 13411 . 



B. Spatial link integral in the strong coupling limit 

In the finite temperature (T) treatment, we obtain the effec- 
tive action of quarks {\, X) and temporal link variable (Uq) by 
integrating out the spatial link variables (Uj). We shall evalu- 
ate the spatial partition function, 



H. EFFECTIVE ACTION 

A. Lattice QCD action 

We start from the lattice QCD action and the partition func- 
tion with one species of staggered fermion \ with a quark 
mass m . Gluons are represented by the temporal link (U ) 
and spatial link (Uj,j = 1, 2, • • • , d) variables, 
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Here the trace of the plaquette Up is defined as, 
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In this action, M x denotes the mesonic composite, M x = 
XxXx, and we have defined two other mesonic composites, 
V ± , which contain the temporal link variables. Sums over 
color indices are assumed. Quark chemical potential on the 
lattice (//) is introduced as a weight of the temporal hopping 
in the exponential form [62], and the staggered phase fac- 
tor Tfj X = (— l) x °^ l^jf-i in the spatial action is related to 
the Dirac's 7 matrices lf25l l32tl . By using a 75 -related fac- 
tor e x = (— l) 20-1 *~ Xd , a stagg ered chiral transformation is 
given as Xx ^ ^ "Xx Il25ll72ll73[l . The lattice kinetic action 

(r s) 

S F ' is invariant under this chiral transformation in the chiral 
limit too — > 0. 

Throughout the paper, we consider the color SU(iV c = 3) 
case in 3 + 1 dimensions (d — 3). Temporal and spatial lat- 
tice sizes are N T and L, respectively. While T = l/N T takes 
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and integrate out the spatial link variables (Uj). In the strong 
coupling region (g ^> 1), we can treat the plaquette action 
term (Sg oc l/g 2 ) through the expansion in the power series 
of l/g 2 (strong coupling expansion). 

In the strong coupling limit (SCL), we can omit Sg, and 
the spatial partition function is decomposed into that on each 
link, 



4cl = J m =11 / dOi exP [-*,>] 



(8) 



We can carry out the Uj integral on each link ( j, x) by utiliz- 
ing the one-link SU(A r c ) group integral formulae [77], 



(9) 

c—eCbd—f j (10) 



dU UabUL = —SadSbc , 

N c 

J dU U ab U cd ■ ■ -U e f = £a 



and other higher order integral formulae. The spatial part of 
the hadronic effective action density is obtained as 02611 . 

J dUjexp[-s jjX ] = exp[-s ( f x ) ] , (11) 



„( eff ) - 
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+ A j , x {B x B x+ ~ j + (-l) N <(h.c.)) 



(12) 



where B x = e a (x°X X c )x/N c \ represent baryonic 
composites. Coefficients (A n , Aj )X ) are summarized in Ta- 
ble [J for N c — 3. The spatial partition function in SCL is 
obtained as, 



7 («) 

J SCL 



exp[~Y, s j,x] 



(13) 



J-x 



The sum over spatial directions ^ . in Eq. ( [T3l would give 
rise to a factor d due to the spatial isotropy. Suppose that the 
action ^ . x Sj x stays finite at large d, the quark field (x, X) 



4 



TABLE I: The coefficients of the hadronic composites in the effective 
action at the strong coupling limit. Detailed explanation to calculate 
these coefficients are found in Ref. 1731. 



Coefficients 


Values (N c = 3) 


Ai 


-l/(4iV c ) 


A 2 


-(iV c 2 • (JV e - 2)! - 7V c !)/(32 • N 2 C ■ N c \) 


A, 


-(2 ■ 7V C ! - iV c 3 ■ (N c - 2)!)/(128 ■ 7V C 4 ■ AT C !) 


Aj,x 





x + OO 



O 

V+ V: 




x + j 



M M - V~V+ .. 



FIG. 1: Effective action terms in the strong coupling limit and 1/g 2 
corrections. Open circles, Filled circles, and arrows show x, X, 
U v , respectively. 



should scale as d 1 ' 4 . As a result, the mesonic hopping term 
J2j M x M x+ j remains finite 0(l/d°), while higher power 

terms of quarks are found to be suppressed as 0{\/\fd) for 
N c > 3. This is called the systematic 1/d expansion, which 
is proposed first in the application to the Ising model |p78ll - A 
spin exchanging term J^j S x S x+ j is assumed to be finite at 
large d, and the mesonic hopping . M x M x+ j could be ana- 
logue of that l3(ill . In the leading order of the 1/d expansion, 
the SCL spatial partition function becomes, 



4cl = exp — M x M x+] + 0(1/ Vd) 



(14) 



In the third diagram of Fig. Q] we display the leading order 
diagram of the 1/d expansion. 



C. Strong coupling and cluster expansion 



by using the cumulant expansion 07911 . 

00 ( -i\n 
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where (•• -) c is called a cumulant, and corresponds to the 
correlation part in the connected diagram contributions, e.g. 
(Si ) c = (Si) - (S G ) 2 . We find that the effective action 
from plaquettes is expressed in terms of cumulants as, 



AS off =-log(e- SG ) = -J2 
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(18) 



The n-th term in rhs is proportional to l/g 2n , and we can iden- 
tify n = 1 term as the next-to-leading order (NLO) effective 
action, and n — 2 term as the next-to-next-to-leading order 
(NNLO) effective action. 

The above identification of the effective action and the 
strong coupling order is consistent with the cluster expansion. 
In the first line of Eq. ( fT71 ), average of Sq is decomposed into 
cumulants as, 
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(S'g) c , 
(19) 



where the sum runs over all partitions satisfying J2 a na ~ 
n. The plaquette action Sg is proportional to a large volume 
factor ^ x ~ Vol. and a small coupling factor 1/g 2 , hence it 
is necessary to count both of them. The first term in Eq. ( fl9] l 
is estimated as, 



(so): 



0(1/^", Vol 



(20) 



In comparison, other terms have smaller powers in volume. 
For example, the cumulant of the n-th power operator is pro- 
portional to (Vol.) 1 , 



In order to evaluate the plaquette contribution Sg, it is use- 
ful to define an expectation value, 



(O) 
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^SCL 



VUjO[Uj]e- b " , 



(15) 



which has a normalization property (1) = 1. The full spatial 
partition function Eq. (0 can be expressed as, 



As) 



(16) 



It is well known that the expectation value of the exponential 
form operator with a small factor {i.e. 1 / g 2 ) can be evaluated 
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{a; i}£conn. 
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2n 



)(l/ 9 2n ,Vol. 1 ), (21) 



where "{x^} G conn." represents that the sum over = 
1, • • • , n} is restricted to connected diagrams, and such a sum- 
mation is C(Vol. ). In a fixed order of 1/g 2 ™, (Sg) gives 
the leading order contribution in terms of the volume in the 
thermodynamical limit, Vol. — > oo. Resumming all lead- 
ing order contributions in volume of the connected diagrams, 
we obtain the exponential form shown in the rhs of Eq. dT8T >. 
This resummation corresponds to the so-called cluster expan- 
sion 117911 . and is consistent with the strong coupling expansion 
of the effective action presented in Eq.[T8l 
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D. NLO Effective action 

The NLO contribution to the effective action Eq. (TT~8T > is 
found to be, 

ASnlo = (S G ) C = -4 E ( u r + U p)c ■ ( 22 ) 
y P 

As long as we consider NLO terms, the cumulant (• • • ) c is 
equivalent to the expectation value (•••), 



(U P ) C = (U P ) = (1/41) / m Up 



e 



(23) 



We decompose the spatial kinetic action Sp to the plaquette 
related and non-related part, 



s P = 



E 



(24) 
(25) 



In the above sum Ylp> f° ur links, (j, (k, x + j), (j, x + 
k), (k,x) are included in a spatial plaquette P = (jk,x) as 
shown in the fifth diagram of Fig. [TJ The link integral in 
Eq. d23b is also decomposed as follows, 



(U P ) C = 
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s (rff) 

e j - x 



(j,x)£P 



dU p Up exp [— s p] 



Y\ e s *>v J dU ktV exp[-sfe !y ] 



(k,y)(P 



(26) 



The second line shows the plaquette non-related part, and 
is found to be unity from the definition of shown in 
Eq. ( fTTT i. The prefactor in the first line corresponds to the 
normalization factor in the plaquette related part. This factor 
is also found to be unity in the leading order of the 1 jd expan- 
sion, since contains four quarks, and is proportional to 
d^ 1 at large d. 

The Up integral part in Eq. (|2"oT l contains at most only four 
links, and we can perform the link integrals by using the group 
integral formulae Eq. ( TTOb again. For a temporal plaquette, 

P = (JO, x), we find EHll, 



(27) 



where the first term in Eq. ( T27l i corresponds to the leading 
order in the 1/d expansion and shown in the fourth diagram of 
Fig-Hi Note that the temporal link variable Uq remains in the 
non-local color singlet composites V^, and will be integrated 
out later. For a spatial plaq uette, P — (jk, x), Up integral in 
Eq. (EU) is found to be (mB, 



where MM MM term is the leading order contribution and is 
illustrated in the fifth diagram of Fig.Q] 

Substituting Eqs. ( T2"7l i and (l28T > in Eq. d22"i i. we obtain the 
NLO effective action from plaquettes, ASnlo- The plaque- 
tte sum J^p m ASnlo leads to J2j ~ d ar, d Tljk ~ d 2 for 
temporal and spatial plaquettes, respectively. Since the quark 
fields (XiX) scales as d -1 / 4 , the composites in Eq. d2Tb and 
C8] ) scale as V~V + ~ d" 1 and MMMM ~ cT 2 . Putting 
all together, V~V + and MMMM give O(l/d ) contribu- 
tions in ASnlo- 

In the following, we consider the leading (SCL) and the 
NLO in the strong coupling expansion, and the leading order 
of the 1/d expansion. The effective action is found to be, 



S e ff 
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where AS< T ' and AS^ come from the temporal and spatial 
plaquettes shown in Eq. d27b and d28l ) including their hermit 

(t) 

conjugates, respectively. The considered contributions S F , 
MM, and AS^ T ' S ' are summarized in Fig.Q] 



m. EFFECTIVE POTENTIAL 

The effective action derived in the previous section still 
contains quark fields (x>X) an d temporal link variables Uo. 
In this section, we obtain the effective potential !F e s by inte- 
grating out these variables in the mean field approximation, 

J V[U , X, x] e- S °« = J V[$] e -^L-^[*] 

« e-^x-^WI . (34) 

I stationary v y 

In this step, several auxiliary fields ($) including the chiral 
condensate a are introduced on the same footing, and the NLO 
effects are self-consistently evaluated. 



A. Effective potential in the strong coupling limit 



(U jk , x ) c 



16iV 4 x x+3 x+ts x+k+j 



O d~ 5/2 



(28) 



Before discussing the NLO effects, we briefly summarize 
the procedure to obtain the effective potential in SCL. The ef- 
fective action Sscl contains the chirally invariant four-fermi 
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term MM. We apply the so-called Hubbard-Stratonovich 
(HS) transformation. The four fermi term MM is reduced to 
bilinear forms in x ar, d X by performing the Gaussian trans- 
formation with a auxiliary field a, 
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where 6 CT = d/2N c , and the matrix 14j 
^k-j »)/^^ represents the meson hopping. In the last line, 
cr is assumed to be a constant, which is determined by the 
stationary condition dJ 7 ^ 1 " /da — 0. Under this condi- 
tion, the auxiliary field a is found to be the chiral condensate 
a = - J2x( M x)/( N r Ld )- Thus the non-linear term MM 
is converted to the quark mass term b a axx, where the finite 
chiral condensate a spontaneously breaks the chiral symmetry 
and generates the quark mass dynamically. 

Now the total effective action reduces to a bilinear form of 
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m q M x 



N T L d —a 2 
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(36) 



where m q = niQ + b a a represents the constituent quark mass, 
and the inverse propagator of quarks is given as, 

= S -f(e»U , x S x+()>y 



' ,,,) +m q S xy , (37) 



We take account of finite T effects by imposing periodic 
and anti-periodic boundary conditions on link variables and 
quarks, respectively. We take the static and diagonalized 
gauge (called the Polyakov gauge) for temporal link variables 
with respect for the periodicity [34], 



C/ (r,x) = diagte^/^, • • • , e f01f '^ Ifr ) 



(38) 



The corresponding Haar measure is given in the form of the 
Van der Monde determinant, 



dUn = 
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ni 
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(39) 



where the delta function reflects the SU(N C ) property, i.e. 
the baryonic effect in the temporal direction. Owing to the 



static property of the auxiliary field a and the temporal link 
variable in the Polyakov gauge, the partition function Z q = 

f Uo e~* G x is completely factorized in terms of the fre- 
quency modes. Hence the quark path integral can be done in 
each mode independently, and leads to the simple product in 
the frequency By utilizing the Matsubara method (see 

for example the appendices in Refs. H4Ql 14411 ). we obtain the 
partition function as, 



n[ 



dtfox e^/ T TT \i + e -(E q -»)/T+ie 



1 + e -( E i+v)/T+i8 



(40) 



where N T 1 is identified as temperature T, and E(m q (a)) = 
sinh -1 [m g (cr)] corresponds to the quark excitation energy. 
Substituting Eq. d39b for Eq. d40l ), the remnant Uq integral can 
be carried out in a straightforward manner Il34ll (the explicit 
procedure is summarized in the appendix in Ref. iPolP . The 
resultant effective potential ll34l Ha] is a function of the chiral 
condensate a, temperature T and quark chemical potential fi, 



F^ h (a-T,fi) = ^a 2 + V q (m q (a);T,ri , (41) 



V q (m q ;T,fx) = -Tlog 



XN c (m q ) + 2 cosh 



N c n 



T 



X Nc (m q ) = 



smh[(N c +l)E q (m q )/T] 
smh[E q (m q )/T] 



(42) 
(43) 



The same result is also obtained by another method l36Tl . by 
utilizing recursion formulae. The phase diagram is obtained 
by performing the minimum search of the effective potential 
F^£ L , and its structure has been investigated in Refs. l37l[38l 

39, 4a mm. 



B. Extended Hubbard-Stratonovich Transformation 

We shall now evaluate NLO correction terms AS^'^ in the 
effective action Eq. d29b in the mean field approximation. The 
temporal plaquette term AS^ T ' is composed of the product of 
different composites V x and V~ The standard HS transfor- 

mation shown in Eq. (135b cannot be applied for such a term. 
Hence we apply here a recently developed method named ex- 
tended Hubbard-Stratonovich (EHS) transformation [45, 46]. 
Let us consider to evaluate a quantity e aAB , where [A, B) and 
a represent arbitrary composite fields and a positive constant, 
respectively. We can represent e aAB in the form of Gaussian 
integral over two auxiliary fields (ip, cf>), 

e aAB = I dLpd( j )e - a {( v ~(A+B)/2) 2 + (4,-i(A-B)/2) 2 }+ a AB 

dtpdcf.e-^-^+^+'t' 2 -^- 3 ^} . (44) 
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The integral over the new fields (p, </>) is approximated by the 
saddle point value, tp = (A + B) /2 and <f> = i {A - B) /2. 
Specifically in the case where both (A) and (B) are real, 
which applies to the later discussion, the stationary value of 
(\> becomes pure imaginary. Thus we replace <fi — ► iaj and 
require the stationary condition for the real value of to, 



,a4B 



-a{<p 2 -{A+B) v -iu 2 + {A-B)u} 



stationary 



■(45) 



In the case of A = B, Eq. (l45l l reduces to the standard HS 
transformation. We find that e a is invariant under the scale 
transformation, A — > XA and B — ► \~ 1 B. In our previous 
work B4411 . a similar invariance exists but is broken after the 
saddle point approximation. As a result, a careful treatment 
is necessary in order to determine the explicit value of the pa- 
rameter. In the present derivation, the scale invariance is kept 
in rhs of Eq. (l45l l. since the combinations ip — u> — (A) and 
p + u) = (B) transform in the same way as A and B, respec- 
tively. This means that the effective potential is independent 
of the choice of A. 

Now we apply EHS to NLO terms. For the spatial 
plaquette action terms, AS^ S \ we substitute (/3 s /d(d 
l),M x M t 



x+j 



m x+% m x+%+ .) in E q- 63' and obtain ' 



AS^ 



1% 



d(d-l) 



E [rf 



x,0<k<j 



(cp s - u s )M x M x+] - (<p s + u s )M x+k M x _ 



-k+ji 



N T L d -(3 s p 2 s 



E M - M : 
x,j>0 



x+j 



(46) 



In the last line, we have assumed that the auxiliary fields take 
constant and isotropic values. Under this constant auxiliary 
field assumption, ui s effects disappear and the sum X)o<fc<j 
leads to a factor d(d — l)/2 for the p 2 s term. As shown in 
the last line, the coupling terms of ip s and M are rearranged 
to the same form as the meson hopping term in the SCL 
effective action by using the translational invariance. This 
M x M x+ jip s term can be absorbed into the meson hopping 
effects M x M x+ j in the SCL, 



b a + 2(3 s p 



2d 



M T M, 



x+j 



(47) 



TABLE II: The auxiliary fields and their stationary values. In the 
stationary value of tp T , (po — N c — Z x fh q + (3 t uJt. 



Aux. Fields 


Mean Fields 


Stationary Values 


a 


<-M> 


-(l/Z x )(dVjdm q ) 


Ps 


{MM} 


a 2 


(p T 


-{(V + -V-)/2) 


2vV(l + \/l+4/3r<A>) 




-{(V++V-)/2) 


—dV q /dfi = p q 



For the temporal plaquette action AS^ T \ we substitute 

(a, A, B) = (/3 T /4d,-V; + (/i),y" and obtain, 

X-\-J 



AS(T) ^ E [^+[^ + (//)-r- ; , /( , 



x,j>0 



x+j 



VT 



(50) 



In the last line, we again assume that the auxiliary fields p T 
and oj t are constant and isotropic, then ^ . leads to a factor d. 
We combine V terms in Eq. d50l ) with those in the SCL tem- 

(t) + 

poral action S F , and the coefficients of V are found to be- 
come, Z T /2, where Z± = 1 + f3 T (p T ±uj T ). We rewrite these 

coefficients as Z± = Z x exp(±<5/Li). Thus the SCL terms Sp" 1 
is modified by the temporal NLO effects as, 



sP = \Y j [z.v+{h)-z+v+{ 1 j)] 

X 



2 



(51) 



Thus the spatial NLO contributions lead to a shift of the coef- 
ficient for the meson hopping effects, which can be evaluated 
by introducing the chiral condensate a via Eq. (l35l l. The coef- 
ficient modification is cared by replacing S a in Eq. d35l l with, 



where, 



S a = N T L d "-^a 2 



(48) 



^Z- 



zl- (52) 



where b a — dj (2N C ) + 2fl s <p s . The constituent quark mass is 
found to be modified as, 

m' q — m Q + b a a . (49) 



In this way, temporal NLO contributions are expressed as the 
quark wave function renormalization factor Z x and the dy- 
namical shift of chemical potential 5fi. In Table [TTJ we sum- 
marize introduced auxiliary fields. 
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C. Effective potential 

Now the effective action reduces to a bilinear form in terms 
of the quark fields (x, x)> 



S cff =Z x ^2xxG X y(m q ,jl)xy 



N T L d 



"a _2 . Ps 2 , Pt I 2 2\ 



(53) 



where G~y(m q , p) is given in Eq. d37l i with modifications 
(m q , p) — > (rrig, /2). We note that the constituent quark mass 
is modified again due to the quark wave function renormaliza- 



tion factor Z 



A' 



to 



(54) 



Equation d6TT l indicates the stationary value of ui T is nothing 
but the quark number density p q . The stationary conditions 
are summarized in Table [ill 

The auxiliary fields ip TjS are found to be explicit functions 
of a and ui T via stationary conditions, while uj t becomes a 
(T,p) dependent implicit function, uj t = p q (cr,uj T ; p,T). 
Hence we need a self-consistent treatment in the minimum 
search of T c s (a, lo t ) in order to determine vacua. This is a 
consequence of the multi-order parameter (a, lo t ) treatment, 
and a new feature compared with the previous works IT361 l37l 

si. 

The auxiliary field uj t may be interpreted as a repulsive vec- 
tor field for quarks. In relativistic mean field (RMF) models 
of nuclei B80I1 . the isoscalar- vector field uj contributes to the 
energy density as, 



e v 



m^u /2 + g ujN p B (p B )u> ■ 



(62) 



The remnant integrals f _ T , e Zx ^ x.G{m q ,jx)x can eva i_ 
uated in the same manner as in SCL. The effective potential 
(free energy density) is obtained as a function of the auxil- 
iary fields $ = (er, ip T>s ,u) T ), the temperature T and the quark 
chemical potential p, 



^ cff ($; T, p) = .F aux ($) + V,(m,($); T, p) 
b a a 2 (3 s ip 2 s (3 T 



. rT / 2 2 

2 



(55) 

N c \ogZ x , 
(56) 



where V q (rh q (&);T, p) has the same functional form as that 
in SCL Eq. d42l except for modifications m q — > TO g ($) and 
p — > /i. The additional term — iV c logZ x , which has no 
counterpart in the SCL, appears from the quark wave func- 
tion renormalization factor Z x through the fermion determi- 
nant contribution, — log[det(Z x G -1 )]. 

We have introduced four kinds of auxiliary fields <i> = 
(<r,<p s ,ip T ,u) T ), and it may contain some redundant degrees 
of freedom. This can be cared by considering stationary con- 
ditions shown in Eq. d45b . 



cir 



9$ 



9$ 



dV q dm 
dfh q 9$ 



, dVq dp_ = Q 



dp 9$ 



(57) 



Note that V q depends on the auxiliary fields via the two dy- 
namical variables m q and p. Substituting a for $ in Eq. d57l ), 
we obtain the relation, 



1 dV q 

Z x drh q 



(58) 



By utilizing this result, the stationary condition for ip s leads to 
ip a = a 2 . Substituting tp T and lu t for $, we obtain a coupled 
equation for ip T and uj t , whose solution is found to be, 



2tpo 



1 + VI + Ww> ' 
<p = N c - Z x m q + /3 t uj 2 



dp 



dp 



(59) 
(60) 
(61) 



where u> is the temporal component of the omega meson field 
ui v . The negative coefficient of u> 2 results in the repulsive po- 
tential for nucleons, +g UJ N^-', and the coupling with the baryon 
density p B leads to the shift of p B as, 



E + g^u: - p B = E - (p B - g u uj) = E - p E 



(63) 



The saddle point constraint gives w a p B . Most of these char- 
acters apply to the auxiliary field uj t introduced in the present 
work. For example, the u> T contribution to the effective action 
in Eq. (T50b is rewritten as —(3 T ut T 2 /2 + (3 T oj T p q , and the sta- 
tionary condition is oj t = p q . When we replace quarks with 
baryons and introduce an appropriate scaling factor for uj t , 
the above two points are consistent with the properties of oj 
vector field in RMF. 



IV. PHASE DIAGRAM EVOLUTION 

In the previous section, we have derived an analytic expres- 
sion of the effective potential JF cff , which contains effects of 
the next-to-leading order (NLO) of the l/g 2 expansion. In 
this section, we investigate the phase diagram evolution with 
the finite coupling effects j3 = 2N C / g 2 based on the effective 
potential T c g. By developing a self-consistent treatment of 
two order parameters, (a, u T ), we study the j3 dependence of 
the critical temperature, critical chemical potential, the critical 
point and the phase diagram. We also discuss partially chiral 
restored (PCR) matter. 



A. Self -consistent treatment in vacuum search 

Solving the stationary condition of JF off with respect to a 
and lu t , 



dTcs dTcs 



da 



dcu T 







(64) 



corresponds to searching for a saddle point of !F e g in the 
(a, u> r ) plane. Since the quark number density p q is an in- 
creasing function of p which is a decreasing function of w T , 



9 



N c =3, p=2N c /g z =4.5 

*eff 








FIG. 2: (Color online) The effective potential as a function of a and 
w T at (7» = (r c , M= o/2,0),(0.2, A t c , T =o.2),(0,^ c ,T=o) in the 
lattice unit. The solid line represent the set of points which satisfy 
the stationary condition Eq. d61| >. 



the stationary condition for uj t , 



p q , has a single solu- 



tion, uj t = lu t (c, T, /i), for a given value of a at finite 
T. The coefficient of oj t 2 is negative in T c s, hence the so- 
lution gives a maximum of T c s for a given a. Thus the sta- 
tionary point in (er, u> T ) is the saddle point of F e s, at which 
.Feff is convex downward and upward in a and ui T direc- 
tions (d 2 f r e s/da 2 > and d 2 J r e g/du> T 2 < 0), respectively. 
Generally, we may have several solutions of Eq. d64"l >, among 
which the lowest T c s dominates the partition function. 

In Fig. [2] we show T c « as a function of (a, lo t ). Solid lines 
show the solution of the stationary condition for uj t , uj t = 
oj T stat ■ , and filled circles show the saddle points. At \i = 0, 
T c s becomes an even function of u> T , and the stationary value 
of bj T is always zero as shown in the upper panel of Fig. [2] 
At finite fi, we have to solve the coupled equations do*4b self- 
consistently. The solution uj T stat - is a smooth function of a 



-2 





T/T =0, 0.2,..., 1.0, 1.2 (jll=0) 



|U,/ju, =0, 0.2,..., 1.0, 1.2 (T=0) 




|3=2N c /g 2 =4.5 



0.5 1 1.5 2 2.5 3 

a 

FIG. 3: The effective potential as a function of a on the T-axis (upper 
panel, fi = 0) and on the /i-axis (lower panel, T = 0) with /3 = 
2N c /g 2 = 4.5 in the lattice unit. The filled circles represent the 
minimum points. 



at finite T as shown in the middle panel of Fig. [2] In the 
case of T = and finite /i, dT c s/du} T is discontinuous at 
fl = E q . This discontinuity comes from the functional form 



of the quark free energy at T = 0, 



V g (m g ,/i,T = 0) 




< (65) 



The ridge found in the lower panel of Fig.|2]corresponds to the 
line Jl = E q , where T c $ takes a maximum value for a given a. 
The stationary value w T stat - at finite T approaches this ridge 
in the limit T — > 0. Thus the stationary condition for w r is 
found to be equivalent to searching for lu t which maximizes 
J-e.fi for each a also at T = 0. 

The effective potential as a function of a for given (T, fx) 
is defined as, J r e ff(o') = •? r eff(o', w = w stat (tr)), whose mini- 
mum point corresponds to the equilibrium. In Fig. [3] we show 
J-"eff(cr) on the T-axis {p, = 0) and on the /i-axis (T = 0) at 
/3 = 4.5, as an example. The chiral phase transitions in these 
cases are found to be the second and first order, respectively, 
as in the case of SCL. In the following subsections, we discuss 
the nature of these phase transitions. 
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2.5 
2 



e 

II 



II 

H 

£ 0.5 
ii 

06 





(3 C (T=1/N X , |i=0)(MC) 
T C (SCL, |a=0)(MC) 




M-c,T=</ T cn=0 



N c =3, NLO 

2 3 ,4 5 
P=2N c /g 2 



FIG. 4: Critical temperature, chemical potential, and critical cou- 
pling. In the upper panel, solid and dashed curves show the NLO 
results of T c , M =o and /i c ,T=o, respectively. Solid squares show the 
results of_u Ci T=o in the previous work by Bilic, Demeterfi and Pe- 
tersson 13711 . We also show the MC results of the critical coupling 
(J3 C ) for given N r = 1/T = 2, 4 and 8 at fj, = (filled triangles) and 
the critical temperature in the strong coupling limit (open triangles). 
From the left I5ll l8lll (the SCL result of T c with monomer-dimer- 
polymer simulations ). 18211 (J3 C = 3.67(2), N T — 2, with a quark 
mass m = 0.025), fmW c = 3.81(2), N T = 2, m = 0.05), (H 
(/3 C = 4.90(3), 8 3 x 4 lattice, extrapolated to m = 0), E3] (J3 C = 
5.037(3), 16 3 x 4 lattice, m = 0.05), 0] (J3 C = 5.040(2), 6 3 x 4 
lattice, mo = 0.05), and 1 85] (N T = 8, extrapolated to mo = 0). In 
the lower panel, we show the ratio R = /x c ,t=o /2c, m =q. 



B. Critical temperature at zero chemical potential 

Along the T axis (// = 0), the quark number density p 9 = 
uj T is always zero, then the effective potential under the oj t sta- 
tionary condition is simply given as !F c s(a) — !F c ff(a,u! T — 
0). In the upper panel of Fig. [3] we show the effective po- 
tential at several temperatures (T/T c = 0, 0.2, . . . , 1.0, 1.2) 
at n = for = 4.5. The effective potential has one lo- 
cal minimum in the region a > 0. As T becomes large, the 
minimum point of the effective potential smoothly decreases 
to zero from a finite value. We find that the phase transi- 
tion along the T axis is the second order as in the case of 

scl gill 13. 

In Fig. 2J we show the critical temperature at zero chem- 
ical potential T c fl= o as a function of 0. We find that T c is 
suppressed as becomes large. This decrease would be a 
natural consequence of finite coupling, since hadrons are less 
bound than in SCL. In the present treatment, the decrease 



of T c is caused by the wave function renormalization factor 
Z x in Eq. d54l ), which has a similar effect to the temporal 
lattice spacing modification. The second order phase tran- 
sition temperature is obtained from the condition C2 = 0, 
where T c s — J2n C n a n /n\. The effective potential !F e s{a) 
and auxiliary fields ^ = ((p TiS ,ui T ) are even functions of 
a in the chiral limit /da \ a ^o = 0, and the first deriva- 
tive of auxiliary fields are zero from the stationary conditions 
dTcfi/d^ = 0. By using these, we find C2 is given as, 



C 2 



d_ 
da~ 



E 



da d^J 



da 2 



bl 9 2 V 



da" 2 

N c {N c + l){N { 



Z£ dm 2 



cr=0 



2) 



Z 2 . 3T(N C + 1 + 2cosh(7V c /VT)) 
From the condition C2 = at \i = 0, we find, 



T, 



T, 



(SCL) 



c,fi— 



Z l 



1 d(N c + l)(N c + 2) 
Z 2 6(N C + 3) 



(66) 



(67) 



(SCL) 

where T c represents the second order phase transition 
temperature at fi = in the strong coupling limit. It should 
be noted that at // = and a = 0, Z x does not depend on 
the auxiliary field. As shown in Eq. d67| i, the critical tem- 
perature decreases due to the wave function renormalization 
factor, Z x > 1 at a = 0. This mainly originates from the sup- 



pression of the constituent quark mass fh q 



m'JZ x . 



In this 



way, the decrease of the critical temperature is understood as 
the rhg modification effects caused by the plaquettes. The T c 
values here are consistent with those in Ref. |39]. 

In Fig. [4] we also show the results of the critical coupling in 
Monte-Carlo simulations with N r = 2 B, 4 |[Toll76ll84ll and 
N T = 8 18511 temporal lattice sizes. These results corresponds 
to T c = 0.5, 0.25 and 0.125. Results with N T = 2 and 
N T = 4 are those with mo = 0.05, 0.025 or O(extrapolated), 
and chiral extrapolated results are shown for N T = 8. Com- 
pared with the results of the critical coupling f3 c = 5.097 
in the quenched calculation with N T = 2 18311 . j3 c is signif- 
icantly smaller with finite masses, B c — 3.81 and 3.67 for 
mo = 0.05 and mo = 0.025 18211 . The monomer-dimer- 
polymer simulations on anisotropic lattice 151 
critical temperature T c = "f 2 /N r = 1.401(2) 
is the critical anisotropy in the chiral limit. The decrease of T c 
in NLO at finite is not enough to explain these MC results, 
and higher order effects such as the next-to-next-to-leading 
order (NNLO) and Polyakov loop effects would be necessary. 



C. Chiral transition at finite density 

At finite /1, the quark number density p q = uj t is gener- 
ally finite and depends on a. We search for uj t which max- 
imizes jF e ff(cr, u>r) for a given a, and we substitute the solu- 
tion, oj r = u> T sta,t -(a), in the effective potential. In the lower 
panel of Fig. [3] we show the effective potential at several 



8 111 give the 
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chemical potentials (/i//i c = 0, 0.2, . . . , 1.0, 1.2) at T = 
for j3 — 4.5. The effective potential has one local minimum 
for fi smaller than the second order critical chemical poten- 
tial, (i < /if, 2 y^ , and two local minima appear in the larger /i 
region. For (3 — 4.5, the vacuum jumps from the NG phase 
(cr ~ cr vac ) to the Wigner phase (a = 0) at critical chemical 
potential, fi = \i c , and this transition is the first order. 

The first order chiral transition at finite /i necessarily in- 
volve the density gap. In the case of (3 — 4.5 and T = 
shown in the lower panel of Fig. 12 the effective potentials at 
two points in (cr, w T ) plane become equal at fx = \i c , and the 
first order phase transition takes place. In the Wigner phase, 
the quark mass is small (zero in the chiral limit), then the 
quark number density is high. At high densities, the chemical 
potential effects are reduced as fi — > fi = fi — 5fi as discussed 
in the previous section. 

In the upper panel of Fig. |4j we show the critical chemi- 
cal potential at T = 0, (J, c ,t=o, as a function of (3. In the 
region of f3 < 6, the phase transition at T = is the first 
order, as in the case of SCL results. We find that the first or- 
der critical chemical potential fJ-^La is not largely modified 



(SCL, 1st) 



0.55. For 



from the strong coupling limit value /i 
example, we find /i lst) ~ 0.58(0.60) at = 3.0(4.5). This 
small modification is understood as follows: In the low tem- 
perature region, the first order phase transition is described in 
terms of the competition between the quark chemical potential 
and the constituent quark mass. Since the temporal plaquette 
suppresses both, the relative relations between them are not 
largely changed. Hence fii 1 ^ ~ /xi SCL ' lst ' ) follows. Results 
by Bilic et al. JHHEHl] are also shown in Fig.gJ Our results 
are qualitatively consistent with their results. 

We can now discuss the critical value ratio R = 
I^c,t=o/T c .^=o, which characterizes the shape of the phase 
diagram. In the lower panel of Fig. [4] we show this ratio as a 
function of (3. As already discussed, T C)/J= o rapidly decreases 
as (3 increases, while the finite coupling effects give rise to 
only small modifications of fi c ,T=o- As a result, the ratio R 
significantly increases with j3 as shown in the lower panel of 
Fig. [4] The ratio R becomes close to 1 at (3 = 6, and much 
larger than the SCL results R SCL ~ 0.3 - 0.45 gUSH. 
The lattice MC results indicate that the critical end point may 
locate in the region fi/T > 1.0, which suggests R > 1.0. 
Based on the recent MC results (T c = 170 - 200 MeV) and 
a naive estimate N c /j, c > Mm (Mm is the nucleon mass), the 
expected ratio in the real world would be R = 1.5 ~ 3. Thus 
the finite coupling effects are found to increase the ratio R and 
make it closer to the empirical value. 



D. Partially chiral restored matter 

One of the characteristic features of the present treatment is 
that the second order critical chemical potential /ii 2nd ^ is finite 
even at T = 0. The second order critical chemical potential 
^(2nd) ^ a gj yen rp j s b tamec j by so i v i n g the condition C% = 
at finite chemical potential. By using Eq. d66l ). the shifted 
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FIG. 5: (Color online) The phase diagram for the several value of 
[3 = 2N c /g 2 in the lattice unit. The solid and dashed lines represent 
the first and second order transition lines, respectively. The actual 
transition is described by the thick dashed and solid lines. 



second order critical chemical potential is obtained as, 



-(2nd) = JL arccosn 
N r 



(N C + 3)(T^ CL) /Z^-T) 



2T 



(68) 



In rhs, a = is assumed. This equation is an implicit equa- 
tion, which should be solved with the condition lo t = p q si- 
multaneously; Z x in rhs is a function of u> T , which is a func- 
tion of jl. The second order critical chemical potential differs 

from /2c 2nd ' 1 by <5/i, 



M ( 2n d) = A (2nd) +Sfl= -(2nd) + lf)g JzJZ. 



(69) 



0, jlc nd ' ) becomes zero. In SCL, we do not have the 



At T 

second term in Eq. (|69l and /i£ 2nd ^ approaches zero at small 
T. In NLO, the second term in Eq. d69l is finite at finite (3. 

As a result, there is a possibility that /ic 2nd ' ) overtakes ni lst \ 
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-1.80 

T=0.09, 0.10, 0.15, T CEP 
T CEP =0.1548, ^ c (lst) (T) 
-1.82 - p = 2N c /g 2 =6 

* -1.83 




FIG. 6: Effective potential on the first phase transition boundary at 
(3 = 6. Solid lines show the results at T = 0.12, 0.13, 0.14, 0.15 
and Tcep, where Tqep = 0.1548, and dashed lines connect the 
coexisting equilibrium points on the boundary. 



which leads to the appearance of the partially chiral restored 
matter 00. 

We numerically find that /i£ 2nd ' can overtake pi lst ^ in the 
case (3 > 4.5. We show the phase diagrams with (3 — 
3.0, 4.5 and 6.0 in comparison with the SCL phase diagram 
in Fig. |5] The tri-critical point (TCP) starts to deviate from 
the second order phase transition boundary at j3 ~ 4.5, and 
becomes a critical end point (CEP) at larger (3 even in the chi- 
ral limit. When the CEP exists off the second order phase 
transition boundary, we have the temperature region, where 
the second order critical chemical potential is larger than the 



first order one, pc (T) > p [ c lsz, (T). In this temperature re- 
gion, T c ff at p = /ii lst ' should have two local minima in the 

a > region as shown in Fig. [6] At p = p{ lbt \ FeS at two 
local minima are equal, and the local minimum with smaller 
a cannot be at a = 0, since the curvature of !F e g is negative at 
around a = (i.e. C2 < 0) in the chemical potential region 

. (2nd) 
P < Pc ■ 



(1st), 



> 



In the temperature region where the condition pc 
pi lst ^ is satisfied, we have three chemical potential regions, 

p < p[ lst \ p[ lst ^ < p < pi ad \ and p > /«£ 2n . The vac- 
uum is in the NG phase in the first region, where the chiral 
condensate is large enough. In the third region, the chiral con- 
densate is completely zero, and it is in the Wigner phase. In 
the second region, the chiral symmetry is weakly but sponta- 
neously broken, and a partially chiral restored (PCR) matter is 
realized 11451 14611 . It is interesting to investigate the transitions 
among them. As p increases, the a jumps from the NG local 
minimum to the PCR local minimum with a > 0, as we can 
guess from the jF e g behavior in Fig. [6] For larger p, the chiral 
condensate in PCR matter decreases, and the Wigner phase 

(2nd) 



fist) 



(a = 0) is realized at p = p c 

The appearance of PCR matter, or equivalently, p y c""' < 
p < /ii 2nd ^ re gion may stem from the multi-order parameter 
treatment 11451 14611 . To clarify this point, we examine several 
truncation schemes, where the effective potential T e ff(a,p q ) 



systematically reduces to that with a single order parameter a. 
And we would check the disappearance of the PCR matter. 

The first treatment is the same as that we have discussed 
in previous subsections, and abbreviated as NLO-A. In the 
second treatment (NLO-B), 0{l/g A ) contributions in Z x and 
p are truncated as, 



^(NLO-B) 
X 

~(NLO-B) 



=1 



=/' 



(70) 
(71) 



In this treatment, we find that tp T and lo t couple to quarks 
separately through rh q and p, respectively. In the third pre- 
scription (NLO-C), we further truncate 0(l/g 4 ) terms in rh q 
and in log Z x . 

^(NLO-C) =({(ff(T + m())(1 _ + 2j3slps(X ^ (?2) 

A^ aux = -iV c logZ x 

» - N c f3 T <p T (NLO - C) . (73) 

It is also possible to expand V q with respect to Sfi = \i — jl 
(NLO-D), 

V( NLO - D )(m g ;^T)~V g (m 9 ; M! T)-/3 TWT ^ . (74) 

These truncation schemes are summarized in Table HTT1 

Stationary conditions in NLO-B, C and D are solved in a 
similar way to NLO(NLO-A). In NLO-B and C, uj t is still an 
implicit function of a, uj t — p q (<7, uj t ; fi, T), and the multi- 
order parameter property is still kept. In NLO-D, lo t is explic- 
itly obtained as a function of a, 



(nlo-d) = dV q (m q (a); fx,T) 
dp 



(75) 



where the rhs does not contain co T . The dynamics is described 
by a single order parameter a. In this meaning, the NLO-D 
gives a similar formulation to those in the previous works i37l 
lillsi]. 

In the upper panel of Fig. [7] we show the phase diagrams 
in NLO-A, B and C in the large p region. In NLO-A and B, 
the maximum temperature of the first order phase boundary 
decreases, and the critical point deviates from the second or- 
der phase transition boundary at j3 ~ 4.5 and 3.0 in NLO-A 
and NLO-B, respectively. In NLO-C, the second order critical 
chemical potential /zi 2nd ' ) at T — overtakes the first order 
one at f3 ~ 3.5. Between the first and second order phase 
boundaries, we find PCR matter. In NLO-D, PCR matter does 
not appear in any region of (T, p, f3). 

In the middle and lower panels of Fig. [7] we show the com- 
parison of p q and a in the present treatments NLO-A, B, C and 
D ll45ll46ll . The gradual increase of the quark number density 
p q after the first order transition is a common feature of the 
multi-order parameter treatments 11451 14611 . This means that p q 
is high, but still smaller than the maximum density iV c in PCR 
matter at small T. At low temperatures, we can investigate 
the appearance of the PCR matter more intuitively. The quark 
number density p q = — dJ r e g/dp is evaluated as, 



Pg_ _ 2smh[N c p/T] 

V ~ X Nr + 2cosh[N c p/T] T^ 1 



,.AV 



(76) 
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TABLE III: The truncation schemes in NLO-A, B, C and D. In NLO-C and D, 0{l/g 4 ) terms in m q are truncated to be m, 

(b a (j + m )(l - P T ip T ) + 2f3 s ifi 3 a. 



(NLO-C, D) 

i 



5fi 



A^ aux 



NLO-A 



lop 



-N c log ^/Z+Z- 



V q (m q ,p,T) 



NLO-B 



1 + (3 T ifr 



-iV c log(l + I3 t <Pt) 



V q (m q ,p,T) 



NLO-C 



/3 T U) T 



^(NLO-C) 



-N c f3 T ip T 



V q (rh q ,p,T) 



NLO-D 



■ (NLO-D) 



BV 

V q (m q ,p,T)~f3 T m T ^ 



where x = exp 
isfied at small T 

E, 



~{E q - p)/T]. When E q > p is sat- 
we obtain x — > and p q — > 0, while 
q < p leads to x — > oo and p 9 — > N c . Medium density 
< p q < N c can appear only in the case where the en- 
ergy and chemical potential balances, E q = p, and x stays 
finite at T = 0. Since p is a decreasing function of u> T , we 
may have a medium density solution of Eq. ( fTol i in the region 
p(tj,u) T — N c ) < E q (<T,ui T ) < p. Specifically in NLO-B and 
C, E q = /i is found to be equivalent to the density condition 
p q = (p, — E q )/f3 T , which can take the a medium value. In 
the large region, this medium density matter can emerge in 
equilibrium and corresponds to the PCR matter as indicated 
in Fig. Q Also in NLO-A, PCR matter appears in a similar 
mechanism at finite T. Thus the multi-order parameter treat- 
ment is essential to obtain the PCR matter at low T, and we 
observe the two chiral transitions as p increases. 

Now we have found following common properties as long 
as the quark number density is treated as the order parame- 
ter in addition to the chiral condensate |{45il46ll . (I) The par- 
tially chiral restored (PCR) matter can appear in the large (3 
region, (II) PCR sits next to the hadronic Nambu-Goldstone 
(NG) phase in the larger p direction, (III) the quark number 
density is high as 0{N C ) in PCR, (IV) in PCR matter, the 
effective chemical potential is adjusted to the quark excita- 
tion energy, and (V) the second order chiral transition to the 
Wigner phase follows after NG^PCR transition. All these 
properties would be the essence of the quarkyonic matter and 
transition proposed in Ref. M86I1 . In the previous work, the 
quark-driven Polyakov loop evaluated in SC-LQCD is shown 
to be small as 0(1/N C ) [36], and it would not grow much 
at low temperatures. This feature is also consistent with the 
proposed property of the quarkyonic matter. 

The quarkyonic matter is originally defined as the confined 
high density matter at large N c [86], and recently investigated 
by using the PNJL model l87ll . In order to discuss the decon- 
finement dynamics, the higher order of 1/g 2 expansion would 
be essential, and a subject to be studied in future. 



E. Phase diagram evolution 

We shall now discuss the phase diagram evolution with (3. 
In Fig. [8] we show dependence of the phase diagram on (3 in 
NLO(NLO-A). As (3 increases, the second order phase bound- 
ary is compressed in the temporal direction according to the 
decrease of 7^ 2nd ) , Phase transition of cold (T = 0) dense 
matter is calculated to be the first order, and the critical chem- 
ical potential does not move much. 

Thick lines in Fig. [8] show the coupling dependence of the 
first order phase transition boundaries. As more clearly seen 
in Fig. [5] the slope of the first order phase transition bound- 
ary becomes negative at finite coupling, and it is natural from 
the Clausius-Clapeyron relation. We expect that the entropy 
and quark number density is higher in the Wigner phase, 



As 



,(W)_„(NG) 



> and Ap q 



4 W) - 



■p q NG ^ > 0, where 



s (w,NG) an( j p( W ' NG ) denote the entropy and quark number 
density, respectively, in the Wigner and the NG phases. With 
this expectation, the slope of the first order phase boundary 
from the Clausius-Clapeyron relation, dp/dT = —As/Ap q , 
should be negative. This improvement from SCL may be re- 
lated to the gradual increase of the quark number density in 
the PCR matter and in the Wigner phase. In SCL, the quark 
number density jumps to an almost saturated value, p q ~ A^ c , 

at low T at p, = p!£ ■ In this case, the lattice sites are al- 
most filled by quarks, and the entropy density will be very 
small in the Wigner phase. This density saturation is a lattice 
artifact, and is expected to be weaken at finite (3. With NLO 
effects, the vector field lo t suppresses the sudden increase of 
p q , and the quark number density gradually increases after the 
first order phase transition. We have discussed this feature in 
the PCR matter, and it also applies to the quark matter in the 
Wigner phase. 

We also find that the slope dp/dT is always negative in 
whole T — p plane at finite (3. This point is different from the 
previous NLO works j39ll . 

The end point of the first order phase transition boundary 
is the critical point, which is either the tri-critical point (TCP) 
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FIG. 7: In the upper panel, solid and dashed curves show the first 
and second order phase transition boundaries, and dots show the crit- 
ical end point. In the middle and lower panels, solid, dotted, dashed, 
and dot-dashed curves show the results of NLO-A, B, C and D, re- 
spectively, and dots and open squares show the first order transition 
points. 



FIG. 8: The phase diagram evolution with the finite coupling effect 
[3 — 2N c /g 2 in the lattice unit. Thin line surface shows the bound- 
aries between the chiral broken and restored phases, and thick line 
surface shows the first order boundaries. 
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or the critical end point (CEP). As we already discussed in 
the previous subsection, the TCP at small (3 deviates from the 
second order phase transition line at (3 > 4.5, and becomes 
the CEP. The temperature of this critical point gradually de- 
creases, while the chemical potential stays in a narrow range 
as (3 increases. We show the evolution of the critical point 
with j3 in Fig. [9] The decrease of the critical point tempera- 
ture, Tap, is consistent with the results in the NJL and PNJL 
models [89]. In these works, it is demonstrated that Tqp de- 
creases as we adopt a larger vector coupling relative to the 
scalar coupling. In the present work, [3 T ui T is regarded as the 
vector potential for quarks, and it grows as j3 increases. 

In the continuum limit, one species staggered QCD would 
become the four flavor QCD with degenerate masses 17511 . 
where the chiral transition is expected to be the first order due 
to anomaly contributions lf88h . The present behavior of the 
critical point shows that the NLO SC-LQCD does not contain 
anomaly effects. 



FIG. 9: Coupling dependence of the critical point. Filled circles 
show the critical points, and solid and dotted lines show the first and 
second order boundaries for ft — and 6. 



V. CONCLUDING REMARKS 

We have investigated the chiral phase transition in the 
strong coupling lattice QCD at finite temperature (T) and 
chemical potential (/i) with finite coupling (f3 = 2N c /g 2 ) ef- 
fects. We have derived an analytic expression of the effective 
potential using one species of staggered fermion in the leading 
(strong coupling limit; SCL) and next-to-leading order (NLO) 
of the strong coupling (1/g 2 ) expansion and in the leading 
order of the large dimensional (1/d) expansion. We have fo- 
cused our attention on the phase diagram evolution. 

From the NLO effective action, we have derived the ef- 
fective potential under the mean field approximation based 
on a self-consistent treatment of NLO effects with a recently 
proposed extended Hubbard-Stratonovich (EHS) transforma- 
tion 1 45, 46]. Then the quark number density (p q ) is naturally 



15 



introduced as an order parameter. NLO contributions are ex- 
pressed via the shift of the constituent quark mass, dynamical 
chemical potential and the quark wave function renormaliza- 
tion factor. The NLO effective potential is found to become a 
function of T, /i, 0, the chiral condensate a and quark num- 
ber density p q . Such a formulation has been essential in order 
to investigate the mechanism of the phase diagram evolution 
with 0. The phase diagram has been obtained by performing 
the minimum search of the effective potential in the multi- 
order parameter treatment. 

The effective constituent quark mass rh q is found to be sup- 
pressed as increases. As a result, the critical temperature T c 
decreases and becomes closer to the Monte-Carlo results at 
(j, = oHIMil, while it is still larger than the MC data. 
The effective quark chemical potential // is also suppressed as 
becomes larger. We have found the small modification in 
the critical chemical potential fi c at low T, In this way, the 
ratio R — fJ> c ,T=a/T c ,ij,=o becomes larger and closer to the 
empirical value. The (3 dependences of T C u=o and jU C) T=o are 
consistent with the previous results 1 39] . The first order phase 
boundary is found to satisfy dfi/dT < at finite (3. This be- 
havior is natural from the Clausius-Clapeyron relation, and is 
different from the SCL results and previous results with NLO 
effects [39]. In the phase diagram evolution, the tri-critical 
point is found to move in the lower T direction. This trend 
is consistent with model results [89]. Partially chiral restored 
(PCR) matter is found to appear in the low T and the large fi 
region with (3 > 4.5. We have shown that the multi-order pa- 
rameter (a, u> T ) treatment is essential in describing PCR mat- 
ter, where the effective chemical potential is automatically ad- 
justed to the quark excitation energy. 

We have discussed the NLO results in the region < 6, ex- 
pecting that the strong coupling expansion is convergent even 
in the region = 5 ~ 6. In the pure Yang-Mills theory, the 
character and strong coupling expansions seem to be conver- 
gent in the region of 2N c /g 2 ~ 2N C for color SU(2) JH 



and SU(3) [20]. For color SU(3), the MC simulations in- 
dicate that the critical coupling r at u = seems to be a 
smooth function of T = 1 / N T H H IH 

and reaches = 5.08 for N T — 8 [85]. When we take into 
account the next-to-next-to-leading order (NNLO) contribu- 
tions in SC-LQCD with quarks for color SU(3), T CjM= o and 
Mc,t=o are found to be very similar to those in NLO in the 
region < 6 B90I1 . These observations suggest that the strong 
coupling expansion does not break down in the region < 6. 
It would be necessary to investigate the NNLO effects on the 
critical point and PCR matter in order to examine the present 
results. 

There are several points to be discussed further. When we 
take into account NNLO contributions, the Polyakov loop can 
appear from two plaquettes. Hence it becomes possible to 
investigate the phase transitions with three order parameters, 
er, p q and the Polyakov loop. In addition, the Polyakov loop 
contributions in the NNLO may modify the dependence of 
T c . The higher order of the 1/d expansion is also an important 
subject to be studied. The baryonic contributions are included 
in the sub-leading order of the 1/d expansion, and would be 
essential to solve a challenging problem: nuclear matter on 
the lattice. 
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